DEFECT MODES AND HOMOGENIZATION OF PERIODIC 
SCHRODINGER OPERATORS 



M. A. HOEFERt AND M. I. WEINSTEIN* 

Abstract. We consider the discrete eigenvalues of the operator H E — —A + V(x) + e 2 Q(ex), 
where U(x) is periodic and Q(y) is localized on M. d , d ^ 1. For e > and sufficiently small, discrete 
eigenvalues may bifurcate (emerge) from spectral band edges of the periodic Schrodinger operator, 
Ho — — A x + V(x), into spectral gaps. The nature of the bifurcation depends on the homogenized 
Schrodinger operator La,q = — V y ■ AV y + Q(y). Here, A denotes the inverse effective mass matrix, 
^— " associated with the spectral band edge, which is the site of the bifurcation. 
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1. Introduction and Outline. Sclf-adjoint elliptic partial differential opera- 
tors with periodic coefficients e.g. the Schrodinger operator with a periodic potential, 
tIh . the time-harmonic Helmholtz equation with variable refractive index, and the time- 

harmonic Maxwell equations with variable dielectric and permeability tensors, play 
a central role in wave propagation problems in classical and quantum physics. The 
spectrum of such operators, characterized by Floquct-Bloch theory [2S1 12S1 H2] , con- 
sists of the union of closed intervals (spectral bands). The cigenstates are extended 
(not localized) and form a complete set with respect to which any function in L 2 (M. d ) 
may be represented. 

In many problems in fundamental and applied physics, periodic media are per- 
turbed by spatially localized defects. These may appear as random imperfections 
in a media, e.g. a defect in a crystal, or in engineering applications, they may be 
, introduced deliberately in order to influence wave propagation [4j [17] . Since the es- 

sential spectrum is unchanged by a sufficiently localized and smooth perturbation 
, (Weyl's theorem, [H]), typical localized perturbations will only introduce eigenvalues 

■ in spectral gaps of the spectrum with associated localized defect modes. 

This paper is concerned with a class of localized (defect) perturbations to a peri- 
odic Schrodinger operator of the form: 



H £ = -A x + y(x) +e 2 Q( £ x) 



where V(x.) is periodic on M. d , Q(y) decays as |y| tends to infinity and e is a small 
parameter. 

Our main result, Theorem 13. 11 concerns the perturbed eigenvalue problem 

H s u £ =fi e u e , u s eH 1 (R d ), (1.1) 

for e positive and sufficiently small. See section [3] for hypotheses on the periodic 
potential, V, and the localized perturbation, Q. 

For e sufficiently small, we prove the bifurcation of discrete eigenvalues into the 
spectral gaps, associated with the unperturbed operator, H = —A + V(x). For 
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any given spectral band edge, we give detailed expansions with error estimates for 
the perturbed eigenvalues and corresponding localized cigcnfunctions in terms of the 
cigcnstates of a homogenized Schrodinger operator 



L, Q = - | ± A, A + Q(y). (1.2) 

Here, Aji denotes the inverse effective mass matrix, associated with the particular 
band edge from which the bifurcation occurs; see Theorem 13.11 Ajt, derivable by 
formal multiple scale expansion (see section!?]), is expressible in terms of the band 
edge (Floquct-Bloch) eigenstate. It is proportional to the Hessian matrix D 2 E bit (k*) 
of the band dispersion function, associated with —A + V^(x), evaluated at the band 
edge = E bje (k*). 



Referring to the schematics of figures ITTT1 and IT721 wc discuss our results. 

• Suppose the inverse effective mass matrix, A, is positive definite and assume 
La,q has an eigenvalue, ca.q < 0. This occurs if Q(y) is a "down-defect" 
(sufficiently "deep" in dimensions d ^ 3) as in figure II. 11 a. In this case, 
Theorem 13.11 asserts the existence of an eigenvalue at E* + e 2 eA,q + 0(e 3 ) < 
E*. 

• Now suppose the inverse effective mass matrix, A, is negative definite and 
La,q has an eigenvalue, eA.Q > 0. This occurs if Q(y) is a "up-defect" 
(sufficiently "high" in dimensions d ^ 3) as in figure II. li b. In this case, 
Theorem 13.11 asserts the existence of an eigenvalue at E* + e 2 eA,Q + 0(e 3 ) > 
E* 

• Fig. 11.21 shows a more general band edge bifurcation when La,q has three 
eigenvalues s^ q < &aq < e AQi the largest of which is degenerate with mul- 



tiplicity three. Theorem 13.11 asserts the existence of five ordered eigenvalues 



at E*+E 2 e%+0(s 3 ), j = 1, 2 and E* +e 2 e% +e 3 ^ ( 3 k) +C(e 4 ), k = 1,2,3. 



1.1. Outline of the paper and overview of the proof. Section[2]summarizes 
the required spectral theory for Schrodinger operators with periodic potentials and 
introduces variants of the classical Sobolev space, H s (R d ), which provide a natural 
functional analytic setting. Section [3] contains the hypotheses on V and Q and the 
statement of our main theorem, Theorem l3.il In scction|3]we present a formal multiple 
scale / homogenization expansion in which we systematically construct bifurcating 
eigenstates and eigenvalues to any prescribed order. In section [5] we prove Theorem 
13.11 In particular, we study the equations governing the correction, ty £ to the iV— 
term multiple scale expansion. 

To obtain error bounds of suitably high order in s, we use a Lyapunov- Schmidt 
approach. Specifically, we decompose the error into Floquct-Bloch modes associated 
with energies lying near the spectral band edge, E#, and those lying "far" from E*: 
ty £ = ^neaj. + ^f m -- ^ncar nas the character of a wave-packet, spectrally supported on 
a small interval with endpoint E* . The next step is to solve for VPf^ as a functional 
of the "parameter" ^noa,., with appropriate bounds. Substitution of ^farl^near] 
the near equation implies a closed equation for ^noaj.. With strong motivation from 
the structure of terms in the multiple scale expansion, wc appropriately rescale, solve 
via the implicit function theorem, and estimate ^ear- The approach we take has 
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E* + e 2 e A , Q , sgn(A) = +1 =S> e A ,Q< 



E t + £ 2 e A ,Q, sgn(A) = - 1 e A ,Q > 



E (k) 



P -Ei(k) *£ 2 (k) 



B (k) " Bi(k) E a (k) 



Fig. 1.1. aj Periodic structure with "down defect" and corresponding localized eigenstate for 
the case of positive definite effective mass tensor, b) Periodic structure with "up defect" and corre- 
sponding localized eigenstate for the case of negative definite effective mass tensor. Below are shown 
eigenvalue bifurcations from band edges of the form E b Qi*) = E%. 



0{e s 



Fig. 1.2. Schematic of band edge bifurcations in the case where the inverse effective mass 
matrix, A, is positive definite. The homogenized operator, L A q, is assumed to have two simple 
eigenvalues and one degenerate eigenvalue with multiplicity three. 



been applied in the context of the nonlinear Schrodingcr / Pitaevskii equation in 

[21 im EMCEE] • 

Previous work for linear Schrodinger operators: Bifurcation of eigenval- 
ues from the edge of the continuous spectrum for Schrodingcr operators with small 
decaying potentials, corresponding to weak defects in dimensions one and two for the 
case of a homogeneous medium or vacuum (V = 0), was studied in |30j . Conditions 
ensuring the existence of eigenvalues in the gaps of periodic potentials were obtained 
in [1] and |13U14| . using the Birman-Schwinger (integral equation) formulation of the 
eigenvalue problem. Homogenization theory was applied to obtain eigenvalues in the 
spectral gaps of a class of periodic divergence form elliptic operators, governing lo- 
calized states in high contrast media in |18[ [8] . An elementary variational argument 
in spatial dimensions one and two, yielding general conditions for the existence of 
discrete modes in spectral gaps of periodic potentials, was recently presented in [26] . 
More general, variational methods can be applied to obtain defect modes which are 
obtained as infinite dimensional saddle points of strongly indefinite functionals; see, 
for example, [IT]. 

Our results concern a particular class of weak defects, slowly varying and of small 
amplitude: e 2 Q(ex), which give rise to defect modes in any spatial dimension. We 
note that the one- and two-term truncated multi-scale homogenization expansion of 
defect modes, which we construct, are natural trial functions for a variational proof 



4 



M. A. Hoefcr and M. I. Wcinstein 



of existence of ground states; see the discussion in Appendix Note also that the 
scaling of the perturbing potential, e 2 Q(ex), also arises naturally in solitary standing 
wave ("soliton defect mode") bifurcations from band edges of periodic potentials in the 
nonlinear Schrddinger / Gross-Pitaevskii equation U6f . 

Homogenization theory has been used to study periodic elliptic divergence form 
operators near spectral band edges in [51 [71 [5] . Homogenization results for the time- 
dependent Schrodingcr equation with a scaling, equivalent to the one considered here, 
were obtained by two-scale convergence methods in [3] ; see also [28l [HE]- In [3] the 
contrast between the scaling we use and the semi-classical scaling is discussed. These 
results establish the validity of the homogenized time-dependent Schrodingcr equation 
on certain finite time scales. The results of the present paper focus on a subclass of 
solutions, bound states, which are controlled on infinite time scales. 

Finally, we mention work on effective classical electron motion in solid state 
physics, derived from the Schrodinger equation for an electron in a spatially peri- 
odic Hamiltonian, perturbed by spatially slowly varying electrostatic and magnetic 
potentials [HI [Ml Hi], in a semi-classical limit. 

Acknowledgments: This research was initiated while MAH was an NSF Post- 
doctoral Fellow under DMS-08-03074 in the Department of Applied Physics and Ap- 
plied Mathematics at Columbia University. MIW was supported in part by NSF grant 
DMS-07-07850 and DMS-10-08855. MIW would also like to acknowledge the hospi- 
tality of the Courant Institute of Mathematical Sciences, where he was on sabbatical 
during the preparation of this article. 

1.2. Notation and conventions. We note that we may, without loss of gener- 
ality, restrict to the case where the fundamental period cell is = [0, l] d . Indeed, let 
B denote the fundamental period cell, spanned by the linearly independent vectors 
{ri, . . . , r^} and define the constant matrix Hr 1 to be the matrix whose j th column 
is Tj . Then, under the change of coordinates x t—> z = IZx, 

— V x • V x + V(x) acting on L^ OJ .(B) transforms to 
-V z .aV z +U(z)= - £ ay + V-( Z ) 

i,j = l 1 J 

acting on L^ eI ([0, where 

a = mZ T , V(z) = V (ft _1 z) , x = TZ^z . 



1 . Integrals with unspecified region of integration are assumed to be taken over 
^ i-e- ff =/ Rd /(x)rfx. 

2. For f,geL 2 , the Fourier transform and its inverse are given by: 

J-{/}(k) = /(k) = J e- 2 - k x /(x) dx, (1.3) 
^ x { 5 }(x) S g(k)= J e 2 - x k 3 (k)dk. 
Thus, T J"" 1 = Id. 

3. = [0, l] d is the fundamental period cell, tt* = [-1/2, l/2] d is the dual 
fundamental cell or Brillouin zone 
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4. 1a( x ) is the indicator function of the set A; x(|k| ^ a) = ljkefi* : |k|^a} 

5. The repeated index summation convention is used throughout 

6. Fourier spectral cutoff: 

X (|V| < a)G(x) = (.F-^dkl < a)T) G = J e 2 " x ' k X (|k| < a)G(k)dk 

(1.4) 

7. T and T _1 denote the Gelfand-Bloch transform and its inverse; see section[2j 

8. Bloch spectral cutoff: 

X(|V| < o)G(x) s T- 1 { ^ x(| • I < aS bb J% {G} (-)p6(x; •)} (x), 

where 6* is the index of the spectral band under consideration, 

9. H s = H s (R d ) is the Sobolev space of order s 



f\\h - 2 W da f\\v ~ = / (!+ |k|*) 2 |/(k)| 2 dk (1 



5) 



2. Spectral Theory for Periodic Potentials. In this section we summarize 
basic results on the spectral theory of Schrodinger operators with periodic potentials; 
see, for example, [29 j [20 | [12] ■ 



Gelfand-Bloch transform: Given / e L 2 (R d ), we introduce the transform T and 
its inverse as follows 

r{/(-)}(x;k) =/(x;k) = ^ e 27 " z - x /(k + z), (2.1) 
T- 1 {/(x;-)}(x)= [ e 2 ™ k /(x;k)rfk. (2.2) 

One can check that T~ l T = Id. 

Two important properties of the transformation T are Td Xj f = (d Xj + 2mkj) Tf 
and (Te 2 " k 7) (x, k) = e 27rik x 77(x, k). It follows that 

(T$(V)e 2 ™ k 7) (x, k) = e 2 " k:c $(V + 2irik) (Tf) (x, k) (2.3) 
T («(■)/) (x, k) = u(x) (T/) (x, k), if u is periodic. (2.4) 



Floquet-Bloch states: We seek solutions of the eigenvalue equation 

(-A + ^(x))u(x) =Eu(x) (2.5) 

in the form it(x;k) = e 27 " k x p(x; k), k e O* where p(x; k) is periodic in x with 
fundamental period cell Q. p(x; k) then satisfies the periodic elliptic boundary value 
problem: 

(-(V + 2™k) 2 +y(x))p(x;k) = £?(k)p(x;k), x e T d . (2.6) 
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For each k e f2*, the eigenvalue problem (|2.6[) has a discrete set of cigenpairs 
{ P(,(x;k), Eb(k) }b^>o which form a complete orthonormal set in L 2 cr (fl). The spec- 
trum of —A + V(x) in L 2 (R d ) is the union of closed intervals 

spec(-A + V) = |J Eb(k). (2.7) 

b^Q, keSl* 

We will study the bifurcation of eigenvalues from the band edge 

E*=E K (K), k. d e {0, 1/2}, j = l,...,d, (2.8) 

with the associated, real- valued band edge eigenf unction 

W (x) = e 27rik *' x p^(x;k*) e L 2 (n). (2.9) 

For example, the lowest band edge is Eo(0) and the associated eigcnfunction is periodic 
Po( x + e j;0) =po( x ;0), j = l,...,dfor the standard Cartesian basis vectors {ej}!- =1 . 
Remark 2.1. 5mce 

w(x + ej ) = e 27 " fe *^w(x) = sj-w(x), (2.10) 

where Sj = +1 if = and Sj = —1 otherwise, the natural function space to work 
in is L 2 symm (n), i.e. f 6 L 2 ymm {Q) if f e L 2 (ft) and /(x + &,-) = s-,/(x). Without 
Zoss of generality, and for ease of presentation, we focus on the case where Sj = +1, 
j = 1, ...,d so that L 2 ymrn (fl) = L 2 er (fl), the space of square integrable, periodic 
functions. This implies that k* = 0. The more general case in eq. \2. 8\) can be handled 
by taking k — > (k — k*) and interpreting values o/k reflected about the boundary offl*. 
The simplicity of E* and the relation VE 1 ^ (k*) = (see Hypothesis H2 in Sec. 
implies that E(k! + k*) can be extended as an even function of k'j = kj — fcj,* for 

j = i,...,d m- 

We will make repeated use of the following self-adjoint operator 

= -A + V(x) - E* : H 2 cr (fl) -> L 2 (n). (2.11) 



Projections % and Completeness of Floquet Bloch states: Define 

%{f}(k) = <p 6 (-; k), /(•; k)> i2(0) = f M^kj /(x; k) dx. (2.12) 

J Q 

By completeness of the {p&(x; k)}b^ 

/(x;k) = J]T b {f}(k) p b (x;k) 

Furthermore, applying T~ l we have 

f( x )=Vf %{f}(k) u b ( X ;k) dk (2.13) 
b>o Jn * 

= T\ (u b (-,k)J) L 2 {Md) u b (x;k) dk, (2.14) 
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where itf,(x;k) = e 2 ' x p;,(x;k). The second equality follows from an application of 
the Poisson summation formula. 

Sobolev spaces and the Gelfand-Bloch transform: 

Recall the Sobolev space, H s , the space of functions with square- integrable derivatives 
of order < s. Since E o (0) = inf spcc(-A + V), then L = -A + V(x) - E (0) is a 
non- negative operator and H s (WL d ) has the equivalent norm defined by 

\\<f>\\ H s ~ ||(I + L )MU 2 
Introduce the space (see, e.g. [27l 19]) 

X s = L 2 (ft*,Z 2 ' s ), 

with norm 



(2.15) 



fx* 



t (l + H^) S |T 6 M(k)| 2 dk. 

6=0 



(2.16) 



Now note that 



/ e 2-k. y rh{ }(k ) (i + Eb (k) - S (0)) f !%(■, k) dk 

yf \T b mk)\ 2 \i + E b (k) - e (o)\ s dk 
E ( 1 + \ b \ ld Y I \%mk)\ 2 dk 



1^ 



(2.17) 



The second to last line follows from the Weyl asymptotics Ef,(k) ~ [15]. Thus we 
have 

PROPOSITION 2.1. H s (W l ) is isomorphic to X s for s ^ 0. Moreover, there exist 
positive constants C\, C2 such that for all <fi 6 H s (M. d ) 



du\\ HS(m < u\\ x ^c 2 \ 



(2.18) 



3. Main Results. In this section we give a precise formulation of our main 
theorem, Theorem 13. II The following are our assumptions. 

HI Regularity. V e Lf^Sl), Q e H a (R d ) for a > d, E H e C 3 (fi*), and 

^e^ cl .(J];C 3 (l]*)). 
H2 Band edge. E* = £k (k*), where k* is an endpoint of the b 1 ^ band such 
that 

(a) E* is a simple eigenvalue with corresponding cigenfunction 
W (x)=.e 2 " k *- X ^(x;k,) ej ff 2 cr (fi) 
and normalization ||iy||.£,2(fi) = 1. 
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a) 



spcc(L A , Q ), sgn(A) = +l 



b) 



spec(L A , Q ), sgn^) 



-1 



Fig. 3.1. Discrete and continuous spectrum of La,q- a) Positive definite effective mass tensor, 
b) Negative definite effective mass tensor. 



(b) VE b JK)=0- 

(c) The Hessian matrix, 

is sign definite. 

H3 Existence of eigenvalue to homogenized equation. 

Introduce the homogenized operator 

La,q S -V y .4V y + Q(y) = - , £—A jl —+Q{y) (3.2) 

jV ^ 

Set sgn(yl) = +1 if A is positive definite and sgn(^4) = —I if A is negative 
definite. Assume La,q has a simple eigenvalue &a,q with sgn^^^g < 
and corresponding cigenfunction Fa,q{y) g H 2 (M. d ); i.e. 

L a ,qF a ,q = e AtQ F A!Q , / Fl Q (y)dy = l, sga(A) e A , Q < 0; (3.3) 



see figure kTTT a) . 

Remark 3.1. For further details regarding the smoothness properties of Eb^ 
and pi,^ with respect to k ; we refer the reader to \29\ \32}j . It can be verified that 
hypothesis H2 holds in one dimension at all band edges U2f and at the lowest band 
edge in arbitrary dimensions \ 19}/ . Band edges with multiplicity greater than one exist, 
e.g. for the separable potential V(x) = 2j=i Vi{ x j)> 2. 

Theorem 3.1. (1) Positive definite effective mass tensor: Assume hypothe- 
ses H1-H3, with sgn(A) = +1. Then, there exists Eq > such that for all < e < eq, 
(jl.ip has an eigenpair /j, £ , u e (x;fi e ) e H 2 (M. d ). fj, £ lies in the spectral gap of— A+V(x) 
at a distance 0(e 2 ) below the spectral band edge having E* as its left endpoint. 

Moreover, to any order in e, this solution can be approximated by the two-scale 
homogenization expansion, see Eq. (|4.35[) . (|5.1[) . with error estimate 

N 

N 

Me _ - E 2 e A . Q - ^ £ >«| < £N+1 °> ( 3 - 4 ) 

n = 3 

for all N ^ 4 and some constant C > 0, which is independent of e. 

(2) Negative definite effective mass tensor: Assume hypotheses H1-H3, with 

sgn(A) = — 1. Then, the statement of part (1) applies, but now fi £ lies in the spectral 
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gap of —A + V(x) at a distance 0(e 2 ) above the spectral band edge having as its 
right endpoint. 

Theorem 13.11 extends to the case where La,q has multiple and/or degenerate 
eigenvalues with bifurcations from band edges with k* / 0, as discussed in the fol- 
lowing two remarks. 

Remark 3.2. General band edge bifurcations: As discussed in Remark 
\2.1l Theorem \3. 1\ generalizes to band edges where k* ^ satisfying eq. 12. 8\) so that 



™( x ) eL 2 symm (Q). 

Remark 3.3. Multiple simple eigenvalues: Note that if La,q has M (finitely 

many) eigenvalues, &aq, j = 1,...,M of multiplicity one, then Theorem \3.1\ avvlies 
directly. Specifically, there exists Sq > such that for all < e < Sq, there are 

eigenvalue / eigenvector branches e — > f/ii" , \'U e (-;/i£^)J . This behavior is shown in 

Fig. \1.2\ with two simple eigenvalue branches with spacings 0(e 2 ). 

Remark 3.4. Branches emanating from degenerate eigenvalues of La.q- 

In spatial dimensions, d > I, the operator La,q may have degenerate eigenvalues, 
e.g. if there is symmetry in Q(y). Suppose eA,Q has multiplicity M. Then, since Ta,q 
is self-adjoint, eA,Q perturbs, generically, to M distinct branches. Thus, applying the 
method of proof of Theorem \3.1\ each degenerate eigenvalue of La.q of multiplicity M 
gives rise to M branches of eigenpairs of H E . The cluster of M distinct eigenvalues 
of H e are within an interval of size 0(e 3 ) about E% + e 2 eA,Q- The j th eigen-branch 
satisfies the error estimates 

N 

N 

M« - - e 2 e A , Q - 2 j < z N+1 C, (3.5) 

n— 3 

for j = 1,2,..., M , all N ^ 4 and some constant C > 0, which is independent of e. 
This behavior is shown in Fig. \1.2\ where an eigenvalue of multiplicity three bifurcates 
from the band edge. 

4. Homogenization and Multi-scale Expansion. We derive a formal asymp- 
totic expansion for the bound state that bifurcates from the band edge into a gap. 
The results of these calculations will be used as an ansatz in the next section [5] to 
rigorously prove existence and error estimates. 

We assume that it e (x;/i e ) satisfies eq. (jl.ip 

[-A + T/(x) +e 2 Q(ex)] Ue = /i £ it £ , (4.1) 
and expand it in an asymptotic series as follows 

00 00 

u e (x;/i e ) = C/ £ (x,y) = J e"C/ n (x,y), ^ e = E* + J e>„, (4.2) 

n— n—1 

where y = ex is the slow variable. Treating x and y as independent variables, equation 
(|4.1j) then takes the form 

( V x + eV y f + V(x) + e 2 Q(y) 1 f/ £ (x, y) = Me [/ £ (x, y) . (4.3) 
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Wc seek a solution U £ (x, y) which is periodic in the fast variable, x, and localized in 
the slow variable, y. Specifically, we assume t/„(x,y) e Lp 0r (f2; H 2 (M. d )). Inserting 
(|4.2j) and (|4.3|) into Eq. (|4.1[) and equating like powers of e we find 

0(e°): L*U = [-A x + V(x) - E m ]U = 0, (4.4) 
Oie 1 ) : LJJx = 2V X • V y £/ + Mi^o, (4.5) 
0{e 2 ): L*U 2 = 2V X ■ V y *7 x + /xiC/i - [-A y + Q(y) - fx 2 ]U Q , (4.6) 



0(e n ) : L*C/„ = 2V X • V y C/„_! + ^U n - X - [-A y + Q(y) - /i 2 ]C/ n _ 2 (4.7) 

n-l 

+ XI ^i Un -i + VnUo, n^z 3. 

J=3 



Viewed as a system of partial differential equations for functions of the fast vari- 
able x, depending on a parameter y, each equation in this hierarchy is of the form 
L*U = <5(x) where <5(x) has the same symmetry as w(x), the band edge state 
(see (|2.9p ). with period cell fl. To solve these equations, we make repeated use 
of the following two solvability criteria based on the Fredholm alternative applied 
to the self-adjoint operators L* and La.q with ker(L*) = span{u>} c Lp Cr (f2) and 
ker(LA,Q) = spanf-F^g} c L 2 (R d ), respectively: 

Proposition 4.1. Let Q e L 2 er (ft), then L*U = Q has an H 2 er (fl) solution if 
and only if 

<v,gy& (a) =o. (4.8) 



Remark 4.1. If k* # 0, then L 2 er (fl) and H 2 er (fl) are replaced by function 
spaces with the same symmetry as u>(x), L 2 ymm (Q) and H 2 ymm (Vl). See Remark \2.1\ 

Proposition 4.2. Let H e L 2 (R d ), then L a ,qF = H has a solution F e H 2 (R d ) 
if and only if 

<F A , Q ,n) LHRd) =0. (4.9) 



4.1. O(e ) Equation. From H2, there exists a unique, real, bounded eigcnfunc- 
tion w e H 2 (Cl) and a simple eigenvalue E* that satisfy 

L*w = [-A x + V(x) — Et]w = 0, \\w\\ L 2 (n) = l, (4.10) 

so that the general solution to Eq. (|4.4[) has the multiscale representation 

C/ (x,y) =w(x)F (y), (4.11) 

for some Fg(y) e H 2 (M. d ) that will be determined at higher order. 



Defect Modes and Homogenization of Periodic Schrodingcr Operators 



11 



4.2. 0(e 1 ) Equation. Applying Prop. [4~T1 to eq. (|4.5p gives the solvability con- 
dition 

2d yj F (w, d Xj w} L2(n) + AiiF = 0. (4.12) 

Since the integrand in the first term, being the derivative of the symmetric function 
w 2 /2, integrates to zero, 

Mi = 0. (4.13) 

Therefore, the general solution for U\ consists of a homogeneous and particular solu- 
tion 

U^y) =w(x)F 1 (y) +2d y] F (y)L- 1 {d X7 w}(x). (4.14) 

where F\ e H 2 (R d ) is to be determined at higher order. 
Remark 4.2. For d = 1, the general solution is 

(r x dx 1 \ 
-x + J ° ■ (4.15) 

Jo w{x>y J 

4.3. 0(e 2 ) Equation. Inserting the expressions (|4.13|) and (|4.14[) into Eq. (|4.6|) 
yields 

L*U 2 = 2d Vj F l d Xj w - C[F ], (4.16) 
where the linear operator C[G] for G e H 2 (M. d ) is 

£[G](x,y) = - Ad X] L^{d Xl w}^)S yj d yi G{y) 
+ w(x)[-A y + Q(y)-M2]G(y). 

Definition 4.3. Define the operator L AtQ : H 2 (WL d ) L 2 (R d ) by 

L A , Q G(y) = «), C[G](; y)) L2(n) = [-V y • AV y + Q(y) - e A , Q )G(y), (4.18) 

where e A ,Q is the simple eigenvalue associated with the eigenfunction F Al Q(y) in hy- 
pothesis H3 and 

A ji = 5 ] i-^{d X] w 1 L^ 1 {d Xl w}) LHn) . (4.19) 

Proposition 4.4. 

A 3 i = -^d k] d kl E H (K)- (4.20) 

We give the proof in appendix [AJ see also [4] . 
Applying Prop. l4~T1 to Eq. (|4.16|) gives 

(w(-),C[F ](;y)) L 2 (Q) = o L A . Q F = 0, fi 2 = e AQ , (4.21) 

is the effective, homogenized equation for Eq. (jl.ip with the effective mass tensor A. 
We have assumed in H3 the existence of the eigenpair F a ^q e H 2 (R d ) and e A- q e 
^ Thus, -Fb(y) = F A .Q(y). 
The general solution for U 2 consists of a homogeneous and particular solution 

C/ 2 (x,y) = w( X )F 2 (y) + 2d y] F 1 (y)L- 1 {d X] w}( X ) + L-^CiFA^i-^y)}^). (4.22) 
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4.4. 0(e 3 ) Equation. Inserting Eqs. ([4TT5]) . (jHH]), and P~2"2l into equation 
(|4.7[) with n = 3 gives 

L*U 3 = 2d y l ,r, w - C[F X ] - H 3 + ^wF a ,q, (4.23) 

where "H 3 is known 



%(x,y) = -2^.a yj i; 1 K[^, Q ](-,y)}(x) + 

2L- 1 {& 3W }(x)[-A y + Q(y) - e A , Q ]^,F(y). 



(4.24) 



By Prop. [Oj Eq. (|4.23[) is solvable if and only if 

L a ,qFi =-(w(-),n 3 (;y)) L2m +(M 3 F AtQ . (4.25) 
By Prop. 14. 2[ Eq. (|4.25[) has a solution if and only if 

» 3 = (F A , Q (-),(w(°),H 3 (o,-)) L2{n) ) L2 ^. (4.26) 
We can now write Fx in terms of F a .q as 

^i(y) =^! Q {-<Mo),H 3 (o,-)> i2(n) + m FA, Q (-)} (y). (4.27) 

With this choice of F\ , eq. (|4.23|) is solvable and its general solution is 

C/ 3 (x,y) = itf(x)F 3 (y) +2d Vj F 2 (y)L- 1 {d X] w}(*)- 

L^faF^-^) + H 3 (-,y) - » 3 w(-)F A , Q (y)}(x), 



(4.28) 



where ^(y) is to be determined. Note also that ^(y), introduced at 0(e 2 ), is to be 
determined. 

4.5. (e n ) Order Equation. Continuing the expansion to arbitrary n ^ 4 from 
Eq. (|4~7|) we have 



L*[/„ = 2d yj F n _ 1 d mj w - £[F„_ 2 ] - H„ + ^nwF A ^ Q , (4.29) 
where is completely determined by all the lower order solutions Ui , I ^ n — 3 

W«(x >y ) = 2^.5 w i, 1 {/:[F n _ 3 ](-,y) +H„_i(-,y) -M«-i^(-)^Q(y)}(x) 



(4.30) 



i=3 



By Prop. 14.11 eq. ()4.29|) is solvable if and only if 

LA,qF n - 2 = -W').^(-.y)>L 2 (Si) +fi>nFA,Q- (4.31) 
Furthermore, by Prop. 14. 2[ eq. (|4.31[) is solvable if and only if 

f*n = (FA,Q(')><™(°)Mor)) L2(n) ) L , (4-32) 
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(4.34) 



With this choice of F n _2 is given by 

F n _ 2 (y)=L A ^ Q {-(w(o),H n (o,.)y L2{n)+f , n F A . Q (-)](y). (4.33) 

Finally, with this choice of F n _2, U n (x,y) is given by 

£4(x,y) = w(x)F n (y)+2d yj F n _ 1 (y)L- 1 {d Xj w}(x)- 

L^{£[F n _ 2 ](;y) +H n (-,y) - (M n w(-)F AtQ (y)}(x). 
Thus we have: 

PROPOSITION 4.5. The first N equations (|4.4|) . (|4.5|) . (|4.6p . . . f|4.T|) are solvable 
with solutions U n (x.,y) e Hp er (Q; H 2 (M. d )), n = 0,1,..., iV uniquely determined up 
to the two arbitrary slowly varying functions F N (y), F N ^(y) e # 2 (R d ) /or N ^ 1. 
These functions are the slowly varying envelopes of the homogeneous solutions to the 
N — 1 th and N th order equations. Moreover, 

N N 

C/fW) = 2 £ n ^(x,y), A*f } = +e 2 e^, Q + 2 e>„, (4.35) 

n— n— 3 

loii/i i/ie particular choice F n _i = F n = 0, is an approximate solution for the eigen- 
value problem (|4.ip (equivalently (| 1 . 1 1) ) wii/i error formally of order s +1 . 

Remark 4.3. 77ie multi-scale form of the approximate eigenfunction given in 
Prop. \4-5\ is used as a "trial function" in Appendix\B\ to give a "quick" variational 
existence proof for defect modes bifurcating from the lowest band edge. We also show 
that a two term approximation (leading order homogenized solution plus first nontrivial 
correction) yields a better estimate for the energy than the one-term approximation 
(leading order homogenized solution). 

5. Proof of Theorem 13.11 Bounds on u(-,/i) — Ue , fJ- — £te ■ To prove 
Theorem l3.11 we introduce the corrections ^(x) and T e to the approximate solution 
displayed in eq. (|4.35[) through 

u(x;n) = U< N \x,ex.) +e w - 1 * e (x), ^ = ^ +e N+1 T e , (5.1) 

Then the error \l/ e satisfies the equation: 

[£. -£ 2 e A ,Q + £ 2 Q(ex)]<J £ (x) = e 2 R e [V £ , T E ](x), (5.2) 

where 

JV 

R e^e T£] = T e U() + ^ . _ qjj^ + £ ^ n + 1 U N _ n _ e QTJ N 

n=l 

(N-2 \ N 

2 £>n+2 + £ N - 1 ^ + T £ J] 
n=l / n=l 

JV-1 AT 

+ ^ £ " XI Mm + l^V + n-m 



n—1 m—n 

N 



(5.3) 



= T e C/ + 2V X • V y U N - QUn-i + 2 Hn+iU N - n + 0(e) 
= T € wFa,q + U*+ 0(e), 

N 

U* = 2V X ■ V y U N = QUn-x + ^ ti n +iU N - n . 
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The leading order multi-scale approximation Uo(x, ex) = Fa q(£x)iu(x) in the 
ansatz Eq. (|5.1|) with w(x) = e 2 ' Tik *' x p{ )H . (x; k*) the band edge eigenfunction sug- 
gests that the dominant contribution to the frequency content of ^ e will be near 
the band edge E* = £& (k*). Therefore, it is natural to decompose \P e into Bloch 
cigenfunctions 

* £ (x)=V/ P6 (x;k)e 2 ™ k - x T b {^}(k)dk, (5.4) 
b=o Jn * 

with associated energies or frequencies £f,(k) where k varies in the Brillouin zone f2*. 
Moreover, we introduce a spectral localization of ^ £ into frequencies "near" the band 
edge and "far" from the band edge 

* £ (X) = *LarM + *LW = T-^Larlx; ") + T^ff^X; 0> 

*Lar(x;k) =n ncar * e (x) s X (|k-k.| <e')r6,{vl/ £ }(k)^(x;k), 

(5-5) 

§L(x;k) = n far * e (x) = 2 X(|k -k.| > e^^Tb^WPfefrk), 

6=0 

where <5„ jm is the Kronecker delta function and the indicator functions are defined as 

X(|k-k*| < e r ) = l{k 6 n*:lk-k»|< e -}(k), x(|k-k*| e r ) = l{kefi*:lk-k«,|s ee --}(k)- 

(5.6) 

Remark 5.1. For our analysis near the band edge, we will use Taylor expansions 
of various quantities about k = k* . Without loss of generality, we will assume that 
k* = which enables a notationally cleaner presentation. See Remark \2.1\ 

We will use the conventions 



L 2 (0) 



&3j0 



where ^ r ^ ear (k) is a scalar and VE'fj^k) is an infinite vector. This decomposition was 
used in |101 [5J [TB] . The parameter r is assumed to lie in the interval 

re (2/3,1), (5.8) 

the choice of which will be made clear later. 

We now apply the Bloch transform to Eq. (|5.2p . project onto the Bloch modes 
Ph(-; k) and use the properties (|2 .3[) and (|2.4[) to find 

[E b (k) -E,- e 2 e A , Q ] % {* E } (k) + e 2 % {Q(e-)* E (-)} (k) 

(5.9) 

= e 2 r6{i? £ [* e ,T £ ]}(k), 6 = 0,1,.... 
We view this as a coupled system of equations for the near and far frequency compo- 
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near: < 



(5.10) 



far: 



(5.11) 



nents ^(k) and ^(k), kefi* 

f (E b ,(k) - E, - e 2 eA , Q )^ car (k) + 
£ 2 x (|k| < e*)T b . {Q(e-)*U(0}(k) 
= £ 2 X (|k|<0[-r b , {Q( £ -)M>L(-)}(k) 

+ T h , {i^ ear + tff ar ,T £ ]}(k)] 5 

(£fc(k) - J?. - e 2 e^, Q )xt|k| ^ e^ b „ 6 )*Lr, h (k) 
+e 2 x(|k|^ £ ^, b )r fc {Q(e-)*L r (-)}(k) 

= e 2 X (|k| £ ^6.,6)[ - T fc {QH*Sea*(-)} (k) + 
_ r 6 {i? e [^car+*fa r ^ £ ]}(k)], 6 = 0,1,2,.... 

5.1. Lyapunov-Schmidt Reduction. In this section, we derive a functional 
representation of the far frequency components in terms of the near frequency com- 
ponents with an associated estimate. After insertion into the near frequency equation, 
a closed system is obtained. 

We use the implicit function theorem to solve the far frequency equations. To 
this end, we observe the following inequalities due to the definiteness of the matrix 
Aji = a^rdkjdktEbtiO) (see hypothesis H2): 

\E bt QL) - E m - #e A ,Q\ = \Ajtkjh + OQ\l\ 3 )\ 2 Ce* > 0, e 2r sc |k|, keft*, 

\E b (k) -E*- e 2 e A . Q \ > C > 0, \b - 6*| > 1. 

(5.12) 

We now have the following existence result 

Proposition 5.1. There exists e > such that for < e < eq, there is a map- 
ping (^,v,e) -» % r [V,H & far : L 2 (R d ) xR„xR^^ such that $% r [* s near , T e ] 
is the unique solution to the far frequency Eqs. (|5.11[) . The mapping ^j ar is C 1 with 
respect to ip and v. The solution of Eqs. (|5 . 11 1) satisfies the estimate 



|| tart nean 



\\H 2 (M. d ) 



< c 



^ jar\_ ^ near' ^ J 



< Ce 2 - 2r (1 + TT + (1 + e JV - 1 T e )||^ ear || L2(Rd) ) 



(5.13) 



for < r < 1 . 

Proof. Since Eq. (|5.11j) is supported on frequencies away from E* , we can divide 
it by E b (k) — E* — s 2 eA.Q- This suggests studying the equivalent equation G = 
where G has components 



G b (k) = x(|k| ^ e r 6 b ^ b ){Mk) + e 2 [E b (k) - E* - e 2 e A , Q ]- 1 



[% {Q(e-M-)} (k) + % {Q(e-W-)} (k) - T fc {.Rty + </>, w]} (k) 



(5.14) 



6 0. 
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Any function cf> e X 2 satisfying G[<f>, s, ip,v] = is a solution of the far equations 
(|5TTTj) with j/j = 4^ oar , v = T e . G is a continuous map X 2 x K x E -> A" and C 1 
with respect to ^ and u satisfying the estimate 



||G[^e,V,«]|Uo<C[l + i; + (l + e 



JV+l-2r„ 



\\x° + 



Note that 



G[0,0,V>,u] = 0. 



J < oo. 
(5.15) 



(5.16) 



The Proposition follows from the implicit function theorem |23j if we can show that 
D iG[0, 0, ip, v\ is invertible. We have 

(L>iG$,e,ip,v]f) b (k) = X (|k| ^ e r 6 b ,. b )f b (k) 

% {Q(s-)f(-}} (k) - / b (k) (2^- 2 e' w+a + e^v) ( 5 - 17 ) 



+e 2 x (\k\^e r 5 b „ b )- 



E b (k) -E m -e 2 e AX 



Therefore, D^G[0, 0, ip,v] = I is invertible. Note that we use the fact that < r < 1 
to conclude that ]im e ->o £ 2 x(|k| e 1 "^ ,&)/[.£;, (k) — E* — s 2 ca,q] = 0. The implicit 
function theorem implies that there exists £o > and a unique , I'f ar [ , J , near> 6 ^ 
satisfying 



G[*L[*L arl T e ], e ,^ car ,T £ ] 



0, 



(5.18) 



for < e < Sq. 

Equation (|5. 18[) is equivalent to 



(*L\(k) 



T 6 {n far [Q(e-)*L(-) - R e in eai + *I ar , T E ](-)] } (k) 



£ b (k) - E* - e 2 e A ,Q 



(5.19) 



We now demonstrate the inequality in Eq. (|5.13[) . Using (|2.18[) , (|5.11[) and the invert- 
ibility of L* — £ 2 eA,Q to obtain 



I* 



far II /J" 2 ( 



sS C 



4' 



la.r 



A' 2 



e 2 C 



x(\-\> £ r K,b)Tb {Q(eo)*f ar (o) - T £ ](o)} (■) 



2-2r 



C 



(1 + b 2 l d )% {Q(e°)*f» - R £ [fy £ , T £ ](o)} (■ 



.V 2 



1 + 6 2 / d 



.V" 



*S s 2 - 2r C ||Q( £ -)*L(-) " fi e [* £ , T £ ](-)ll i2(R .) 
e 2 - 2 ^ r (1 + £ iv-i T£) (|^ ar( . ) || i2(B<i) + ||^ oar (.)|| i2(R<i) ) + 1 + T £ ] , 

(5.20) 



where the constants C are independent of e. The third inequality results from the 
Weyl eigenvalue asymptotics [15] and the bound (|5.12j) . The last inequality results 
from direct estimation of the error terms (|5.3p . With e small enough so that e 2 ~ 2r C ^ 
1/2, we can subtract the term involving e 2_2r C|| x l , f ar || ff 2( R d) from both sides of the 
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inequality and then divide by 1 — e 2 ~ 2r C(l + e N ~ 1 T £ ) to obtain the desired estimate 
(15331) . □ 

Remark 5.2. Note that we do not obtain smoothness of^fj ar in e. When applying 
the implicit function theorem in the above proof, we did not use any smoothness of the 
map G in e. This is because of the sharp, e dependent cutoff function x(|k| ^ e r ). 

Remark 5.3. The estimate for ^ £ far e H 2 (R d ) and ^ £ near e L 2 (R d ) in (f5T3f can 
also be proved for $f ar 6 H s (R d ) and V £ near e H s ~ 2 (R d ) for s ^ 2. The proof for 
the case s ^ 3 involves application of the operator {I -\- L^Y 1 1 which can be shown 
to be equivalent to the H s ~ 2 norm, to Eq. (|5 ,2|) and necessitates further regularity 
conditions on the functions V(x), Q(ex), and i/ n (x, ex), n = 0, 1, . . . , N. 

5.2. Near Frequency Equation and its Scaling. We now study the near 
frequency equation (|5.10p with the aid of certain Taylor expansions for |k| < s r , 
where we invoke our regularity hypothesis HI 

£*,(k) = E* + ^d kj dkiE h ^{Q)kjki + -dk j d kl d km E bt (k')kjkik m , 

= E* + Ajikjki + -dk j dk l dk m E b ^(k')kjk l k m) (5-21) 

Pb* (x; k) = p b% (x; 0) + d kj Pb* (x; k")fcj 

= w(x) + d k .p bjt (x;k")k j , (5.22) 

for some |k'(k)|, |k"(k)| < e r . Inserting these expansions into the near frequency 
equation (|5.10[) . we have 

[A jl k J k l -e 2 e A ^ £ ncal .(k) 

+ e 2 x (\k\<s r )T b , {Q(e-)^car(-)}(k) =e 2 J RL ar [*Lar,T E ](k), 

where 

^earL^near, T £ ](k) = X (|k] < S r ) [%, {tf [*Lar + *far[*near, ? £ ]} (k) 

-T H {Q( £ -)*L[^ca rI T £ ](-)}(k) 

- ^klkrndu^du^^l^)], (5.24) 

= X (|k| < e r )%* {T £ w(.)F A . Q (e-) + U*(-,e-)} (k) 

+ O x2 [(e + ^+^^I^Larll^], 



where the leading order behavior comes from the definition of R £ in Eq. (|5 .3[) . Recall 
that £/# is defined in Eq. ()5.3|) . The terms proportional to e 3r ~ 2 put a further 
restriction on the exponent r. In order to keep i£J ear order one, we require 

2/3 < r < 1. (5.25) 

In this case, the error term satisfies the estimate 



near I near ; A J 

< C [l + T £ + (e + e 2 ~ 2r + e Sr - 2 )\\n ea rU*) 



(5.26) 
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Equation (|5.23[) can be rewritten suggestively as: 



+ x 



(5.27) 



< e r - L ) {Q(£-)*Lar(-)} (k) = ^ e ar[*aear, T e ](k), 

Thus, we seek a solution of Eqs. (|5.23|) . (|5.24[) in the form 
K cm .(k) = X (|k| < O 1 $ = X < e- 1 ) ^ (7) > ( 5 - 28 ) 

*L ar (x) = r- 1 |x(| ■ I < O 1 * (-) P6s (x; ■) J (x) 



In the next two Lemmata, Lemma 15.21 and 15.31 we express the terms of ()5.27l) , 
which involve the Gelfand-Bloch transform, in terms of the classical Fourier transform 
plus a remainder, estimated to be small in e. 

Lemma 5.2. 

(A) Assume \l/^ eor (k) is given by (|5.28p . Then, 

[ Ajl j. h _ eAiQ ]^ eor ( k ) 

= {( A 3l d yj d m - e A . Q ) X (|V| < e^ 1 ) $} (|) 

76,{Q(e-)*Ur(-)}(k) = p-MQx(|V|< e r ~ 1 )*}^f) + f(k) 

(5.29) 

w/iere 

\\e\\ L * m *S Ce s \\Q\\ H , m m L * m , (5.30) 

with s > d and < r < 1 . 
Proof of Lemma I5.2t Recall the notation Tj = f for the Fourier transform given 
by ()1.3j) . By (|5.28|) since k is localized near we have, Taylor expanding p&^xjk) 
about k = 0, 

*L ar (x)=p 6 ,(x;0)^- 1 |x(H<e' ) ± $ (-) J (x) + &(x) 

= p t ,(x; 0) X (|V y | < e^ 1 ) $(y)| y=ex + £i(x), 
||^i||L 2 (R<i) ^ C£ r |l^lU a (R'')i C = C {\\VkPb^\\L^(n X fi*) d )■ 
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Since T commutes with multiplication by a periodic function (see (|2.4[1 ) and since 
p b# (x; 0) is periodic 



(5.31) 



T{Q(e-)*L ar (-)H*>k) 

= P6,(x;0) T{Q(s-) X (|V e .| < e^" 1 ) S( e .)} (x, k) + T{Q(e.)^(')} ( x > k ) 
By the definition of 7", (|2.1|) . we have 

T{Q(e-) X (|V e .| < e^ 1 ) $(£■)} (x, k) (5.32) 
= ^x {Q(e-) x(|V £ .| <e'- 1 )$(e-) }(k) 

+ 2 1 { X (|°| < e r ) ^ $ Q) j (•) } (k + m) e 2 — 

- T' T ? {Wx(|V| <£•-)*(•) )(i~ 



2 /«(;) is(^H)x(| k + .»-i|<«').^ 



(5.33) 



|m|^l 



To prove (|5.30[) we need to estimate the sum in (|5.33|) for |k| < e r . For such k 
the sum can be estimated as follows: 



m|3:l 



1 



k + m - 1 



X(|k + m-l| < e r ) dl 



< 2 I 



(- 



[m[^l 



I 111 

I e 



2\ 2 



4 dl 

e d 



\m-\\^2e T 



4 <a 



i + 



f dl ||*|[l*{r<«) 



« ce s 2 i m i ' \Q\ii'p*) mwp.*) e s wqwhs^) m L 2 (Rd) , s >d. 

|m|^l 

(5.34) 



A similar calculation shows 
||r{Q(£-)^(-)}lk=(R^fi*) ^Ce s+r \\Q\\ H , 



s>d, < r < 1. (5.35) 



Since T b ,{/}(k) = <p 6 ,(-,k),r{/}(-,k)> L2(n) , we have by Ql, Q, 
JEM]), and (jQ5)) that 

7i.{Q(e-)*L»r(-)}(k) = ^y{Qx(|V|< O*}^) + 5 ( k )' 
II^IU 3 ^") s$ Ce'llQIff.^d) ||$||z,2 (Md) , s>d, 0<r<l, (5.36) 
which is the assertion of Lemma [ 



Lemma 5.3. The right hand side of eq. i5.23\) , defined in eq. \5.21$ , satisfies 

RLaAKear^W (5.37) 



< 



-r-1 



+ ( w 



V e //L 2 (n) 



+ <S(k), 
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where 

||<S|U W ^ C (e + s 2 - 2r + e 3r - 2 \\) mvp*y (5.38) 



Proof. The proof follows in a similar manner as to Lemma l5.2l bv use of eqs. 
and (|5.26l) along with the estimates in Prop. I5TT1 and eq. (|5.24[) . □ 

5.3. Solution of the Near Frequency Equation for Small e. The results 
from the previous section enable us to complete the proof of Theorem 13.11 

Substitution of the expansions of Lemmata 15.21 and 15.31 into (|5.27p and defining 

k 

k — — 5.39) 

£ 

results in 

Proposition 5.4. The near frequency component $ satisfies the equation 

{( Ajtd^dy, - e A , Q ) X (|V| < e^ 1 ) $} («) 
+ Jy{Qx(|V] <e r - 1 )$}(«) =7- y ff[$,T £ , £ ]. (5.40) 

The right hand side has the following form 

FH[^,T^e\ = X (\k\ < O^^eorlxdVI < e^ 1 )^, T s ] + i? c [$]) 

= X(|«| <£ r - 1 )[^F AiQ (K)+(w(-),U*(-,K)) L2 ^ (5.41) 
+ 0{e r + e 2 ~ 2r + e 3r ~ 2 )]. 
We define the following operators Xe > Xe where 

X^xdVyl^" 1 ), Xe^l-Xe = x(Ny\>£ r ~ 1 )- (5.42) 
In physical space we can write (|5.40[) as 

XsL A ,QXs$ = XeH[<i>,r^e]. (5.43) 

where 



JT[$, T £ , e](y) = T e F A , Q (y) + <«,(.), C/ # (-, y)> i2(n) + h[*, T e , e], 



(5.44) 



In order to solve Eq. (|5.43[) . we require a regularization that guarantees the in- 
vertibility of the operator XeL A .QXs- Since zero is an isolated eigenvalue of La,q, 
there is a small disc of radius p = p e about zero, with boundary C p such that for 
e sufficiently small, C p encircles m eigenvalues of XeL Aj QXe, counting multiplicity, 
where m is the multiplicity of zero as an eigenvalue of L A q. 

Introduce the projection onto the spectral subspace associated with eigenvalues 
of XsL A ,QXe, encircled by C p : 

He = [ (XeL A .QXe ~ A/)" 1 dX, (5.45) 

2m Jc„ 
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Note that 

IL =(F AtQ ,-) L2{Rd) , (5.46) 

projects onto the kernel of La,q- 

We now rewrite (|5.43p as the following system for $ and T e : 

XeL AtQX s$ = Xe(I~ Tie)XeH [*, T £ , e] (5.47) 
X£ n e x e F[$,T £ ,£] = 0. (5.48) 

Any solution ($ e ,T £ ) of ([5^47] . (jOSj) is a solution of ([543]) 

We claim that for e small (|5.47|) can be solved for $ = <!> e [T e ] via the equivalent 
nonlocal "integral" equation: 

& = (XsLA^QXe)- 1 (I ~ He) (T'fA.Q + «), [/#(•, y)> i2(Rd) + T £ , e]) . 

(5.49) 

Indeed, the solution may be constructed using the iteration: 

= (x e iA,QX e )- 1 Xe(/-n e ) Xe (T £ F A , Q +< U ;(-),(7 # (-,y)> i2(Rd) +/ l [^,T £ , £ ]), 

^ = ( X£J L A , Q xe)- 1 xe(/ - n £ )(T^ Q + <«,(•), c/ # (-, y)> L2(Rd) ). 

(5.50) 

By use of (|5.45j) and (|5.44j) , we have 

< ||(Xe^QXe)- 1 Xe(/-n e ) Xe (/ l [$ J ,T E ,e]-/ l [$ J _ 1 ,T £ , £ ])|| 

sc r(e) - <Vi II , r(<0 = P" 1 ^ + £ 2 - 2r + £ 3r - 2 ), I < r < 1. 



Therefore, — =SS t^Y \\®i — ®o\\ an d if £ satishes the smallness condition 



t(s) < 1, the sequence {& e Aj^o is Cauchy in H 2 (M. d ). It therefore contains a subse- 



quence, which is convergent to a limit $ £ e H 2 (M. d ). By .ff" 2 (]R d ) continuity of the 
terms in the iteration (|5.50[) . one can pass to the limit in (|5.50[) to obtain a solution 
$1 which satisfies Eq. (|5.49[) . 

This solution $| is a functional of T 6 , and appears in equation ()5.48j) . which we 
view as an equation for T e . We write (|5.48|) in the form 



(5.51) 



g[v,e] = ux £ n e F AiQ (y) + Xe n e (w(-), U#(; y)> i2(n) 
+ XeReh[${v\,v,e] = 0. 

For £ = 0, this equation has the solution g[vo, 0] = with 

vo = - (F A , Q C), <»H ^ # (°, O)^)^ ■ (5-52) 

The Jacobian, D^gfuo, 0] = 1. By the implicit function theorem [23 , for |e| sufficiently 
small there exists a unique solution e i— » T e satisfying <7[Y e ,e] = 0. This completes 
the proof of Theorem 13.11 

Appendix A. Effective Mass Tensor. 

In this appendix we prove Proposition 14.41 relating the Hessian matrix of the 
band dispersion function Ei, (k) to the matrix A resulting from the multiple-scale 
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analysis. In addition, we prove hypothesis H2(b) under certain conditions and the 
positive definitcness of I — A. 

The solutions to the eigenvalue equation (|2.5[) are sought in the form Wb(x;k) = 
e 2 ™ kx p b (x;k), k e fi*. Then, p b ^ and E b ^ satisfy 

4* V* = ("A - 47rik • V + 4ir 2 |k| 2 + V(x) - E bjf (k)) p 6 , (x, k) = 0, x e T d , 

. (A : 1} 

with periodic boundary conditions p bie (x + e^; k) = p b (x; k). Taking the derivative 
of eq. (|A.ip with respect to kj gives 

L ( * ] d kj Pb* (x; k) = (4md Xj - 87r 2 fc, + 8 kj E b ^ (k)) p b:fc (x; k). (A.2) 

Evaluating eq. (IA.2I) at k = k* and using the fact that the kernel of L* k *^ is spanned 
by (x; k*), we arrive at the solvability condition 

dkjEb* (k*) = -4m((d Xj + 2irik* j) p b ^ (•; k*),^ (•; K)) L 2 (n) , (A.3) 

for 7 = 1,2,..., d. When k* = 0, p b ^. (x; k*) is real valued so that eq. (|A.3[) simplifies 
to 

d kj E H (k*)=0, j = l,2,...,d. (A.4) 

For the case d = 1, we also have 

=0, fc* e {0,±l/2}, d = l. (A.5) 

This result follows from properties of the Floquet discriminant A(E) [12]. Briefly, for 
each E, one constructs a 2 x 2 fundamental matrix of solutions M(E) and considers 
the values of E for which M(E) has an eigenvalue 1 or —1 corresponding to a periodic 
or antiperiodic eigenvalue, respectively. This is equivalent to A(E) = +2 where 

trace(M(J5 6 (A;))) = A(E b (k)) = 2cos(27rfc(£; b )). (A.6) 

Differentiating this expression with respect to k and evaluating at k = gives 

^| (fc)) (k) = -47rsin(2^). (A.7) 

Since dA/dE(E) = if and only if E is a double eigenvalue (Theorem 2.3.1, [12]) and 
£7^ is assumed simple (hypothesis H2(a)), eq. (|A.5[) follows. 

The above discussion proves hypothesis H2(b) at the left band edge k* = for 
arbitrary d and both left and right band edges fc* 6 {0, +1/2} when d = 1. It is 
possible for VE b ^ (k*) = in other cases, e.g. separable potentials, and we continue 
the discussion assuming this to be true. 

It follows that 

3 fcj P6*(x;k») =4™(4 k * ) )~ 1 {(a c3 +2™fc* J > b J-;k*)}(x). (A.8) 
Differentiating eq. (|A.2[) with respect to ki and setting k = k* gives 

L^* 1 4 3 -5feiP6*(x;k*) = 

-16tt 2 (^ +27rtfc.,i)(4 k * ) )" 1 {( 3 *j + 2 ^ fc *J>b*(-; k *)}( x ) (Ag) 
-16tt 2 (5 aj +2 7 rifc* J ) (4 k * J ) _1 +2nfc, il ) Wt (-;k t )}(x) 
+ (d^fc^bjk*) -87T 2 (5 ii )p(, 4 (x;k*). 
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Invoking the solvability condition and using \\p bjj . \\L 2 (fi) = 1> CC L- (|A.4[) gives 

+ 16tt 2 ((8 Xl + 2iriK, l )(L i ** ) y 1 {{d Xj + 2irik» d )p b ^ (o; k*)} (■; k,)) 

+ Wn 2 ((8 X] + 2mkt ij )(L < ?* ) )~ 1 {(dx l + 27rifc,,j)Pfc,(o;k,)} (.),P6,(-;k») 



i 2 (n) 
(A.10) 

Integrating by parts and using the fact that (L^. k *^) 1 is self-adjoint, the last two 
terms are equal giving 

-2,2i: 2 ({d X] +2nik*, 3 )p b ^K), (L^)' 1 {(d Xl +2mk^) Pb ^o; K)} (■)) . 

\ / L 2 (Q) 

(A.ll) 

In order to identify eq. (|A. 1 1|) with the final result, eq. (|4.20p in Prop. 14^41 we use the 
definition of w(x) in eq. (|2 .9[) to compute 

d Xj w(x) = e 2 ™ k *' x (d X] + 2mk.j) PbA^M (A.12) 

which is the first term in the inner product of eq. (|A.12[) . In addition, the identity 

L*e 27rlk *' x /(x) = e 27rik *- x Ll k * ) /(x), (A.I3) 

implies 

e 2 ^*-*^)- 1 {(d Xl +27rifc*,0p 6 ,(-;k*)}(x) =L- 1 {d Xl w(-)}(^), (A.14) 

and the result follows. 

Appendix B. Homogenization and Variational Analysis. 

The existence of a bound state for eq. (jl.ljl bifurcating from the lowest band edge 
•Eb(O) can be proved by showing that the Raylcigh quotient 

Pri Jr. {|V^(x)| 2 + [F(x) + g 2 Q(ex) - £ (0)]Kx)| 2 dx} 

£M = /,Mx)Prfx ' (R1) 

is negative for some choice of u e H 1 (M d ) [21]. A natural choice for u is the multi- 
scale expansion in eq. (|4.2j) with e sufficiently small. Furthermore, a higher order, 
two-term trial function gives a better approximation of the energy than the one-term 
trial function. 

Proposition B.l. 

1. Negative energy trial function: with assumptions H1-H3 in section\3\and 
setting E* = Eq(0), the lowest band edge, then there exists Eo > such that 
for allO < e < eo, 

£ [F A , Q (eo)w(o) + 2V F A , Q (eo) ■ L* 1 {Vw} (o)] < 0. (B.2) 

It follows that there exists a ground state. 
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2. Estimates of the ground state energy: if Ljq = — A y +Q(y) has a sim- 
ple eigenvalue eiq < and corresponding eigenfunction FjQ(y) e H 2 (M. d ), 
then there exists Eq > such that for all < e < sq, 

£ [F A , Q (eo)w(o) + 2V F A , Q (eo) ■ L~ l {Vw} (o)] < £ [F I . Q (eo)w(o)] < 0. 

(B.3) 

For the proof of Prop. IB. 11 we will make repeated use of the following averaging 
lemma. 

Lemma B.2. Let p(x) = p(x + z) be periodic with fundamental period cell Q 
and Xizsz d \p*\ < 00 w here p z are the Fourier series coefficients o/p(x). If G e 
L 1 (R d ) n C n (R d ), then 



[ p(x)G(ex)dx -- d \ p(x)dx f G(y)dy 



< Ce n . (B.4) 



Proof. Expand p in the Fourier series p(x) = 2zez d PzC 2,riz ' x . Then 

/ p(x)G(£x)dx = ~ d YpJ G(y)e 2 ™^dy 
jR d e ./r<* 



e 



| f G(y) dy + - d VpJ G(y)e 2 ™^ dy ( B .5) 

JR d £ z ^ JR d 

= ±-f P (x)dx[ G(y)dy + 0(e n ), 
£ Ja Jm d 

where the first line is justified by the assumed absolute convergence of the Fourier 
coefficients and the second line results from integration by parts n times. □ 

First we consider the ansatz u(x) = Fi > q{£x)w(x) for eq. (|B.3[) . A computation 
and several applications of the averaging lemma IB. 21 give 



_ 2 / R . i {|v^Q(y)l 2 + Q(y)I^.Q(y)l 2 }^ , ^ 



J Rd \F LQ (y)\ 2 dy (B.6) 
= e 2 e /iQ +o(e 2 ) < 0, 

for e sufficiently small. 

A similar, more involved computation for the ansatz u(x) = F Ai q(£x)w(x) + 
2VF AlQ (ex) ■ L* 1 {Vw} (x) leads to 

£ [F AtQ (eo)w(o) + 2VF A , Q (eo) - L, 1 {Vw} (°)] 

2 f Rd {VF A , Q (y)-A\7F AQ (y)+Q(y)\F AQ (y)\ 2 }dy 

f R d \F A , Q (y)\ 2 dy +0{£) 1 j 

= e 2 e A , Q +o(e 2 ) < 0, 

for e sufficiently small. For a bifurcation from the lowest band edge, the effective mass 
tensor A is positive definite [19| hence the eigenvalue e A ,Q is negative. 

The proof of Prop. IB. II is completed if we can show that e A ,Q < &i,q- For this, 
we use the following proposition. 

Proposition B.3. The operator I — A = \4(d Xj w, {d Xl w}y L2 ^A is positive 

definite at the lowest band edge (6* = 0, k* =0). 
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Proof. Recall that ^ with one-dimensional L 2 (T d ) kernel spanned by w. Let 
v = («!,..., Vd) G K d be arbitrary. Then 

v ■ {I-A)w = <v ■ Vw,L-\ ■ Vw) L2(n) > E J_ E J v ■ V^|| 2 ^ C||v|| 2 , C> 0, 

(B.8) 

where -Ei(O) — i?* > is the second eigenvalue of L* acting on L 2 (T d ). □ 
Introducing the energy functional 

J K£i .9 2 (y)rfy 

we observe 

, n 7 r i [Vff(y) ■ (A - I)Vg(y)] dy 

by the negative definiteness of A — I from Prop. IB. 31 Using cqs. (|B.6I) . (|B.7[) . and 
(|B~10t we find 

e A,Q = Ja,q[F a ,q] = inf JA,Q[g] ^ Ja,q[Fi,q] < Ji,q[F^q\ = e/, Qj (B.ll) 
and the proof is complete. 
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